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Abstract — This study investigated Human Immunodeficiency 
Virus (HIV) models of Liancheng and Micheal. The model was 
extended by incorporating the treatment term p and a 

1 

polynomial of the form which gives information about 

1 + CfjV 

the current and past states of the virus. The optimality system 
was derived and an optimal control model of drug treatment of 
HIV infection of CD4 + T-cells was investigated. Conditions for 
the optimal control were considered using Pontryagin’s 
maximum principle and solve numerically using Adomian 
Decomposition Method (ADM). Global stability of the 
equilibria, the existence and uniqueness of the solution to the 
problem for the optimal control pair were established. The 
model exhibit two equilibria, disease-free and endemic 
equilibrium. The simulated optimal control pair ( U X ,U , ) 
controls the percentage effect of the chemotherapy on the 
CD4 + T-cells and represents the efficiency of the drug treatment 
in inhibiting viral production and preventing new infections. 
The characterized objective function based on maximizing 
T-cells and minimizing the cost of chemotherapy treatment was 
in agreement with the existing literature. The numerical results 
obtained using ADM was also in complete agreement DTM. The 
result obtained also show that the information on the current 
and past states of the virus would minimize the endemicity of the 

Index Terms — HIV infection, Stability, Basic reproduction 
number, Lyapunov functions, Adomian decomposition method 

I. Introduction 

AIDS is one of the deadliest epidemics in human history. It 
was first identified in 1981 among homosexual men and 
intravenous drug users in New York and California. Shortly 
after its detection in the United States, evidence of AIDS 
epidemics grew among heterosexual men, women, and 
children in sub-Saharan Africa. AIDS quickly developed into 
a worldwide epidemic, affecting virtually every nation. The 
United Nations Program on HIV/AIDS (UNAIDS) estimates 
that the worldwide number of new cases of HIV infection 
peaked in the later 1990s with more than 3 million people 
nearly infected each year. However, some regions of the 
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world, especially Vietnam, Indonesia and other countries in 
Southeast Asia, continued to see an increase in the 
early2000s. In addition, the number of people living with HIV 
or AIDS has continued to rise as a result of new drug 
treatments that lengthen life. 

AIDS is the final stage of a chronic infection with the human 
immunodeficiency virus. There are two types of this virus: 
HIV-1, which is the primary cause of AIDS worldwide, and 
HIV-2, found mostly in West Africa. Inside the body HIV 
enters cells of the immune system, especially white blood 
cells known as T cells. These cells orchestrate a wide variety 
of disease-fighting mechanisms. A particularly vulnerable to 
HIV attack are specialized “helper” T cells known as CD 4 + 
T-cells. When HIV infects a CD 4 + T-cells, it commandeers the 
genetic tools within the cell to manufacture new HIV virus. 
The newly formed HIV virus then leaves the cell, destroying 
the CD 4 + T-cells in the process. 

No existing medical treatment can completely eradicate HIV 
from the body once it has infected human cells. The loss of 
CD 4 + T-cells endangers health because these cells help other 
types of immune cells respond to invading organisms. The 
average healthy person has over 1,000 CD 4 + T-cells per micro 
liter of blood. In a person infected with HIV, the virus steadily 
destroys CD 4 + T-cells over a period of years, diminishing the 
cells’ protective ability and weakening the immune system. 
When the density of CD 4 + T-cells drops to 200 cells per micro 
liter of blood, the infected person becomes vulnerable to 
AIDS-related opportunistic infectious and rare cancers, which 
take advantage of the weakened immune defenses to cause 
disease [9]. A model for the interaction of HIV with CD 4 + T 
cells that considers four populations: uninfected T cells, 
latently infected T cells, actively infected T cells and free 
virus [6]. Effect of AZT on viral growth and T-cell population 
dynamics were considered and numerical bifurcation 
techniques were used to map out the parameter regimes of 
these various behaviours and they showed that when the 
endemic state is stable, it is characterized by a reduced 
number of T cells compared with the uninfected state, thus 
T-cell depletion occurs through the establishment of a new 
steady state. 

An epidemic model of HIV infection of CD 4 + T-cells with 
cure rate and delay were examined [8]. The dynamics showed 
that is completely determined by the basic reproduction 
number Ro<l. If Ro<l, the disease-free equilibrium is 
asymptotically stable and the disease dies out. If Ro>l, a 
unique endemic equilibrium exists and is globally stable in the 
interior of the feasible region. Moreover, they proved the 
effect of that delay on the stability of the equilibria and 
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showed that the introduction of a time delay in the 
virus-to-healthy cells transmission term can destabilize the 
system, and periodic solutions can arise through Hopf 
bifurcation. Model of cell-free viral spread of human 
immunodeficiency virus (HIV) in a well-mixed compartment 
such as the blood stream were studied and discussed the 
existence, stability of the infected steady state and introduced 
a discrete time delay to a CD 4 + T-cell and the emission of viral 
particles on a cellular level [12] .They studied the effect of the 
time delay on the stability of the endemically infected 
equilibrium; criteria were given to ensure that the infected 
equilibrium is asymptotically stable for all delay. 
Mathematical models were considered for the infection of 
human immunodeficiency virus-type 1 (HIV-1) with target 
cells between initial infection, assumed that the infection 
among the cells can be approximated and the classical 
mathematical model with nonlinear infection rate [36]. 

They proved, if R 0 < 1 , then the HIV infection is cleared 
from the T-cells population; if R 0 > 1 , then the HIV infection 
persists. Discrete time delay to the model to describe the time 
between infection of a CD 4 + T cells and the emission of viral 
particles on a cellular level were studied [40]. They examined 
the effect of the time delay on the stability of endemically 
infected equilibrium; criteria are given to ensure that the 
infected equilibrium is asymptotically stable for all delay. 
Also they observed that time delay does not induce instability 
and oscillations in the model. Mathematical model based 
upon the assumption that actively — infected helper T-cells 
come only from the latently-infected T-cell population [37] 
and they studied the time delay that appears in the equations is 
assumed to be the average time that a T-cell remains latently 
infected. 

An optimal control model of drug treatment of HIV infection 
of CD 4 T cells were considered [42, 43] .They showed that 
the optimal controls represent the efficiency of drug treatment 
in inhibiting viral production and preventing new infections. 
Existence for the optimal control pair is established and the 
Pontryagin’s maximum principle is used to characterize these 
optimal controls. The optimality system was derived and 
solved numerically. The results also showed that the optimal 
treatment strategies reduce the viral load and increase the 
uninfected CD 4 T cells count, which improves the quality of 
life of the patient. Series technique for optimal control of HIV 
infection dynamics was reported [41], They obtained the 
minimum amount of medicine required and showed that the 
number of uninfected CD 4 T cells is maximized. Hence, 
minimum amount of medicine ensured that the effects of the 
chemotherapy are minimized. 

n. Model Formulations 

Considered a three dimensional model which consists of 
concentration of susceptible CD4 + T cells (T), concentration 
of infected CD4 + T cells (T*) and free HIV viruses with a 
general non-linear incidence rate. 


— = s-ccT + rT\ 1- 
dt 


T + T 
T~ 


i-pT* (2.1) 


dT 


“— = ^-pr*- P r 

dt 1 + aV 


( 2 . 2 ) 


- = njst* - yv 


(2.3) 


Where T, T and V denote the size of compartmental of 
concentration of the uninfected CD 4 T 
Cells, infected CD 4 T and free virus respectively, s is the 
rate of supply of CD4 + T cells from precursors in the bone 
marrow and thymus, OC is the death rate of uninfected CD4 + T 
cells, k is the rate CD4 + T cells becomes infected by free virus 
while T max is the maximum CD4 + T cells population level. [3 
is the death rate of infected CD4 + T cells, N is the number of 
free virus produced by lysing a CD4 + T cells, y is the death 
rate of free virus, r is the rate of growth for the CD4 + T cells 
population. Qualitative investigation of the system describes 
by Equations (2. 1-2.3) reveals that the long-term behaviou, 
falls into two categories: endemic or dies out. When the 
diseases dies out naturally, the solution asymptotically 
approaches a diseases free equilibrium s Q 
of the form, 

(: r - a)T ^ + ^T m m [(r - a) 2 T mdx + 4rs\ f 
2 r ’ 

The threshold that determines the stability of this equilibrium 

is the R 0 , 




kT 0 Nj3 


(fi+p)r 

When the disease free equilibrium is unstable, there exists an 
endemic equilibrium of the form S x is given as 


e 1 =(7\r ,Vj 


(2.4) 


such that 

f _ f\y + a l Nj3f*)(j3 + p) 
kNj3 

_ - b±Jb 2 -4 ac 
2 a 

y = MEl 

7 


Where 
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P<*1 | 2 f5al p | a\p 2 | ra x P | ra*p \ 

k 2 k 2 k 2 T^k T^k J 


(T r 2 ty 2 p r 2 p ] 

L _M 

r 2Pya x | Aya x p | 2ya ] p 2 \ 

k 2 N 2 p k 2 N 2 p 2 \ 


v k 2 N k 2 N k 2 Np ) 


aa x P aa x p ra x P ra x p ry ryp 

k k k k kT mdx kN kNpT m 

ay ayp ry ryp 

C ~ S ~kN~kNp~~kN + kNp 

Obviously Equation (2.4) will only exist provided > 1 


III. Positivity of solutions 

Theorem 3.1:_Let the initial data T (0) > 0, T*(0) > 0, V (0) > 0, then the solutions T (t), T* (t), 
V ( t ), of HIV free model (2.1) are positive for all t > 0 [45] 

Proof: It is clear from the first equation of model (2.1) that 


dL> _| 

f rT kV 

\t 

dt \ 

^T^l + ay) 

r 

<E>J 

r rT* kV " 


T \ 

K a + T^ + l + a x Vj 

Integrating both sides of (3.2) 


InT > - 

at ~i\ 7- r ( r )- fcy ( r )p 

Take exponential of both sides (3.3) gives 


T{i) > C 2 e L 

Applying the initial condition 


1+aMrU 


T(t)>T( 0)e 




For all r > 0 

Proof: It is clear from the second equation of model (2.2) that 
dT* 
dt 

So that, 


->c e+ p y 

lat, 

T\t) > T* (0) exp |^- £(/? + p)dz^ > 0 
For all r > 0 

Proof: It is clear from the third equation of model (2.3) that 

^->(p+ r y 


So that, 

V(t) > V(0)exp[- £(A + 7>fe] > 0 
For all t > 0 

4.1 Global Stability of Virus Free 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 
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We shall prove the global stability of the virus-free equihbrium by means of Lyapunov function of the system of Equation 
(2. 1-2.3) 

Theorem 4.1: The virus-free equilibrium is globally asymptotically stable if R 0 < 1 . 

Proof: Let us consider the Lyapunov function 


Nj 8 

= -pr-pr* + -§- (Npr* - r v) 

1 + ay NJ 3 


(4.1) 


(4.2) 


1 + ay v r NJ 3 NJ 3 


(4.3) 


1 + ay NJ3 


(4.4) 


( Py{P 

+ r) f 

' NpkVT 


Npil + ay)^) 

[np{i+ 

ay)\ 

+ 

p+p J 

Pr J 

( Mp 

+ r) f 

m T 


Npil + ay)^,) 

[np{ 1+ 

ay)\ 

Jy{P+y) 

P+p) 

Pr ) 




-pr' 


(4.5) 


(4.6) 


< f Py(P + y) [ y K _ 1 + a, y _ Np{\ + ay) r 


\Np(\ + ay){p R ° P + p 


Py 


(4.7) 


<0 

If L =0 then 


f Pr(P+r ) 

r Npk T l + «! 



kPy{P+y) P+p J 

1 Pr J 


(4.8) 


= Oand pT* = 0 


Hence, the largest invariant set included in \L = 0} is reduced to the virus free equilibrium. Thus by Laselle’s invariance 
principle, the VFE is globally asymptotically stable. [44, 46, 47] 


4.2 Global Stability of Endemics Equilibrium 

We shall prove the global stability of endemic equilibrium obtained by means of Lyapunov’s direct method. The Lyapunov 
function constructed of the suitable combinations of composite quadratic and logarithmic functions 

Theorem 4.2 If R 0 >1 then the unique endemic equilibrium S ] of equation (2. 1-2.3) globally asymptotically stable in the 
interior of Q . 

Proof: Define L : {(r,r\v)e Q : T,T*,V > 0 }->9? by 


l(t,t*,v)=^\\t-t |+| r-r |+| v-v 


T*-T* 


P + 2a 

~~k 


(4.9) 
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Then L is C 1 . on the interior of Q , £ x is the global minimum of L on Q and L\T, T , V ) = 0. The time derivative of L 


computed along solution of equation (2.1-2.32) is 


f j3 + 2a V V-V ^dV 

X k ){ V Jdt 

= [(r-r)+(r*-r*)+(y-y)]L-ar + rr(i-^^j-(i-A()/?r t -^ + 


(4.10) 


Using 

s = aT + rT^-X^^ + (l-N)/JT* -yV 
k,vr —(j3 + p)r* =o 


NJ3T* -]V = 0 

To write this, we have 

£’ = [(r - t)+ (r - t‘)+ (v - V)] |- (a + r|(r - r^i - L±U j + {v - v)j + (<* + rfr' - f*)j 

+C8+2a^(r-7-^-^lJ + (r--7-)j + ^±^.(v-v)-(r + (i-v))-r(v-v) (4.12) 

L' = ( a + r|(r-r)^i-^f^j + (v'-v>)j -(r+(i-N))fl(r -fj -^^{v -vj (4.13) 

Hence L 1 is negative. Note that, L 1 = 0 if and only if T = T, T =T and V = V 

Therefore the largest compact invariant set in j(T , T L 1 = 0 } is the singleton } , where £ 1 is the endemic 

equilibrium. LaSalle’s invariant principle then implies that £ 1 is globally asymptotically stable in the interior of Q . [44, 46, 47] 

5.1 Optimal Control Formulation 

We introduced optimal control to the equation (2. 1-2.3) then we have new system of the equation below: 
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dT -s aT l rrfl T + T *] c - 

h kVT , 

1 + pT 

l+a x V 

nt 0 )=r 0 

(5.1) 

dt l ) 

dT* u kVT 

= e 1 BT - pT 

dt 1 + ay 


T\t 0 )=r; 

(5.2) 

dV 

— = e~ u -NBT -yV 
dt 


V(t 0 ) = v 0 

(5.3) 


The model describes the dynamic interactions among healthy CD 4 + Tcells, infected cells and free virus in the organisms and it is 
given by Equation (5. 1-5.3) above. The influence of the antiretroviral drugs in the model Equation (5.1 -5.3) can be simplified 
by applying. Taylor series decomposition as follows: 

e~ u = 1 -u + higher order term. (5.4) 

Then we obtain: 


- = s-aT + rT\ 1 - 


V 1J l + aV 


dt 1 + a,V 


T(t 0 ) = T 0 
T*(t 0 ) = T* 


(5.5) 

(5.6) 


= (l - u 2 )NJ3T* - yV 

V(t 0 ) = V 0 

(5.7) 

I 

1 

s' 1 

II 

1 ^ 

Wit o) = w 0 

(5.8) 


5.2 Optimality system 

Pontryagin’s Maximum principle provides necessary conditions for an optimal control problem. This principle converts 
problem (5. 5-5. 8) into a problem of maximizing and Hamiltonian, H, point wisely with respect to Ui and u 2 : 


H = T(t ) — - [b x u1 + B 2 u 2 2 ] + Y, dJi 


(5.9) 


Where / is the right hand side of the differential equation of i-th state variable. We define the Hamiltoniann function given by, 


H(-) = T(t)- — u 2 - — u 2 2 + Aj s-ctT + rT\ 1- 


-d - Ul )-^-+ P r 

l + a x V 


-pT*y^-u 2 )Npr* -yv)+AXu 1 Npr* -yw) (5.ii) 


Solving the first order condition 


dH dH 

= 0 and = 0 

du t du 2 

(5.12) 

Then we obtain 


f = -^-^ = 0 

du l 1 + a,V 

(5.13) 

— = -B 2 u 2 + (X 4 -y)NpT* = 0 
du , 

(5.14) 


Then the optimality conditions from (5.12 and 5.13) we obtained; 
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kvr 

M. =Z .A., 


" (l + a 1 V)B l 

Npr * .. . . j 

u 2 = — — 

b 2 J 

We obtain optimal expressions U , and u 2 to finally we have, 
kVT 


0 (l + a,V)B 1 




o< / . (^-^)<1 


— a 4 -^> 0 +«,v) 5 , 




(i + «y)s, 

NJ3T 




(A 4 -2 3 ) 




a 4 -^) 


o<®L(a 4 _i 3 )<i 


jspgr* 


1 5, 

However, the optimal conditions is given by 


(/l 4 -^ 3 )>0 


(5.15) 


(5.16) 


(5.17) 


w* = min| 

l,max| 

5 

l 

< 

S 1 

o' 

u 2 = min| 

l,max 

HfH) 


Utilizing equation (5.18) in (5. 5-5. 8) we have the following optimality system 


- = s-aT + rT 1 


+pT 


• It In kV 

1 1 (! + «,’ 


y)B j 




kvT 

(l + «,v) 


dt 

dV 


dW _ 
dt 


, max jo, 
_,max|o, 


5, 


(5.18) 



(5.19) 

] 1 kVT 


■IW)-"-" 

(5.20) 

}}}” 

(5.21) 

\nj3t* -yw 

(5.22) 


Applying (5. 18) in (5. 1 1) and differentiate w.r.t to T, T ,V,W we have; 
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^7 = -l-4| a + r\ |-| l-min<jl, 


4 " " "Up? 


^f = l \ ~K^r~ \ + p\-A 2 (P + p) + ^\ 


y) 

1 b 2 


-^3) 




NP 


= -Aj\ 1 - mini l,max 


x ;„ .v/r ^ 4) |; 


l 4 |l - min jl, max jo. 


NpT\\-ip \\\ kT{\-V) 

(1 + ay)- 


-hr 


8H 

= 

dW 


\np 


(5.23) 


(5.24) 


(5.25) 

(5.26) 


r(0) = T 0 , T\ 0) = r 0 8 ,y(0) = v 0 , w(0) = w 0 , Mtf)=o 

for i = 1,2, 3, 4 

(5.27) 

However, we solve the adjoin equations for problem (5.5-5. 8) 

Using Pontryagin’s Maximum Principle; 

d\ _8H d\ _ -8H cUj _8H dAj _ -8H 

~dt ~ 5T ’ ~dt~ 8T* ~di~^V ~dd~ 8W 

Applying equation (5.28), we have 


(5.28) 



kV 

(5.29) 

1 +ay 

^ = - 1 - ^ [i j + p\ - ^{p + p) + 4,(1 - u 2 )Np) + A 4 u 2 Np 

(5.30) 

dA 3 , ^ ^7^(1- v) , .*r(i-v) 

dt 4 2 (1 +ccyf * 1 (1 +ay) 2 


(5.31) 

^ = -(-A 4 w) = A 4 W 
dt 


(5.32) 

with transversalty conditions 



■si; 

II 

O 

II 


(5.33) 


Theorem 5.1: given optimal controls u i ,u 2 and solution ,V and W 0 f (Corresponding state system of problem 
/ .| . A 2 5 -^3 and A 4 


(5.9), there exists adjoint variable 
by (5.18) 


4 satisfying the Equations (5.33). However, the optimal control is given 


Proof: the adjoint equations and transversality conditions can be obtained by using Pontryagin maximum principle such 
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4 = -^.V', ) = 0; 

The optimal control U 1 and 


4 = -fp.-W»,)=o| 

* = -%' W = °] 

can be solved from the optimality conditions 


du 2 


8H _ 
du x 

That is 

d ^- = -B l u l + {\ + X 1 ) t^- = 0 

OMj 1 + ay 

— = -B lU n + (K + y)NpT % = 0 

SUr. 


(5.34) 


(5.35) 


(5.36) 


By the bounds in u of the controls, it is easy to obtain U 1 and U 2 in the form of Equation (5.10) 

5.3Existence of an Optimal Control Pair 

The existence of the optimal control pair can be obtained using a result by [47]. 

Theorem 5.2: Consider the control problem in (5. 5-5. 8) there exist an optimal control pair 

(«;,«*)£ u such that j(£ tl £) =™ Ui)e J(u v u 2 ) (5.37) 

To proof the theorem in (5.37) an existence result is used and the following properties are checked: 

(1) The set of controls are corresponding to state variables is nonempty 

(2) The control U set is convex and closed 

(3) The right hand side of the state system is bounded by a linear function in the state and control variables. 

(4) The integrand of the objective functional is concave on U. 

(5) There exists constants c l ,c 2 > 0 , and B > 0 such that the integrand L(x, Z, U l ,U 2 ) of the objective functional 
satisfies 

L(x,z,m 1 ,m 2 ) < c 2 -c^l 2 +\u 2 \ 2 ) % (5.38) 

In order to verify these conditions, we use a result by [43]to give existence of solutions of system equation (5 .5-5. 8), which gives 
condition 1 . The control set is convex and closed by definition, which gives condition 2. Since our state system is bilinear in 
ui,u 2 , the right hand sides of system equations (5.5-5. 8), satisfies condition 3, using the boundedness of the solutions. The 
integrand of our objective functional is concave. 

Now to the last condition needed. 

Where c 2 depends on the upper bound on x, and c l > 0 since B X ,B 2 > 0 . Hence, there exists an optimal control pair. 

5.4 Uniqueness of the Optimality System 

Theorem 5.3: The function U ~ triin(max( c, a), b) Ljp SC p,j^ z continuous in c, where a<b are some fixed positive 
constants. 

Proof: consider cl, c2 real number and a,b are fixed positive constants. We will show that the Lipschitz continuity holds in all 
possible cases for max(c,a). Similar arguments hold for min(max( c, a), b) as well. 
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(1) q > a, c 2 > a: |max( q , a) - max( c 2 , a) | = |(q - c 2 )| 

(2) q > a, c 2 < a : |max( q , a) - max( c 2 ,a) | = |(q - a) | < |(q - c 2 1 

(3) q < a, c 2 > a : |max( q , a) - max( c 2 , a)\ = {a - c 2 ) < |(q - c 2 )| 

(4) q <a,c 2 <a\ |max(q,a) -max(c 2 ,a)| = |(a-a)| = 0 < |(q -c 2 )| 


Hence |max( q , a) — max( C 2 , a)| = |(q — C 2 )| and we have Lipschitz continuity of u* in c. 


Theorem 5.4. For sufficiently small final time (tf), bounded solutions to the optimality systems, (5.19-5.22) and (5.29-5.32) are 
unique. 

Proof: With assumption in [49]; suppose (T,T , V, W, A l ,A 2 ,A 3 ,A 5 ) and 

(T,T ,V ,W , A^,/ i 2 ,A 3 ,A 5 )are two different solutions of our optimality system (5.19-5.22) and (5.29-5.32). 

Let 

T = e M p, ,T =e M p*,V =e M q,W = e M x,\ =e~ M w,/ 1 2 = e~ M z,A 3 = e~ M v,A 4 = e~ M y and 
T = e M p, ,T* =e M p*,V = e M q,W = e M x,\ =e~ M w,/ 1 2 = e~ M z,A 3 = e~ M v,A 4 = e~ M y 


Where A>0 is to be chosen. Further we let 


q* = minj 

Jmaxj 

(w-z)ke At } , ] 
a '’(i +a ,v)B t pq \ A \ 

(5.40) 

u* 2 = minj 

|maxj 

p .y\ 

(5.41) 

and 




q* = min< 

|max< 

\ (w - z)ke M 11 

(5.42) 

u* = min< 

max < 

L ( «VW 

(5.43) 

Now we substitute T = e M p into the first ODE of (5.28) and get 


p' + Ap =e M s-ap + rp\ 1 - 


- l-min< max 


xjq, 


(w - z)ke M 

(1 +«,vK 



ke At pq 
1 + a^vq 


P( 5-44) 


Similarly, 

T* = e A, p*;v = e At q, W = e At x, \ - e At w, 1 

(5.45) 

A^ = e^z',^ = e M v,A 4 = e M y \ 


p + A p = 1-mimmax^ q, 


(w - z)ke 


. 3 1 , .1 1 (v-y)Nfie‘ 

q +Aq = 1-mirKmax^q,-^ — — p 


x' + Ax = mirq max 


xjq, 


(l + a x v)B 2 
(v - y)Nj3e' 


-w: 


a,vq 

\N(¥-yq 


~(P+P)p 


(l + ay)B 2 


p \A\ n Pp ~m 


(5.46) 

(5.47) 

(5.48) 
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(v - y)Nj3e M 
25, 



■‘h 


\ke A ‘p - 


— y' + Ay = —yw 


(5.49) 

(5.50) 


From equation (5.29-5.33), we obtain; 

w' + Aw = e' 


- w^a + f j^i 


*-£)} 


%W-z)k e>a pq \ b V\_kpq_ 
(l + ay )B 1 f Mil + cc lV q 


-z' + Az = -vt| 1 + pw-z(P + p)+ l-min| maxj a v V ^ ^ Nfie h p J >,b 2 }- \vNfi + 


"’{f 




yNfi\ 1-mirK max 


Now we subtract the equations for T and T ; T* and T ; V and V ; W and W ; 
A l and A L m ,A 2 and A 2 ;A 3 and A 3 ;A 4 and A 4 

Then multiply each equation by appropriate different of functions and integrate from 0 to tf. 
Next we add all eight integral equations and will use estimates to obtain uniqueness. Using 
Theorem (5.3); we have 

I * _*i 1 I ke h pq ke'pq I 

M, — M, < 7 v — \(W-Z ) (W~Z) d 

1 1 (1 + ay )B l | 1 + a x q 1 + a x q \ 

And |t<2 - m 2 * | - -^-|(v ~ y)Nj3e M p* - (v - y)N/3ke A, p\ 

One case of the estimate was illustrated using on equation tq — u x Hence, we obtain: 


l(p- P) 2 (t f ) + Aj (p - pfdt <aj\ p- p\\p - p\dt + 
2 oo 


(5.51) 


(5.52) 


(5.53) 


(5.54) 

(5.55) 


f|wi - «T|| P - p\dt + p\\p* - p*\p - p\dt < q j[p - pf +| p* - p*\ + 1 q- qf + \x -xf + 

0 0 0 

\w-w\ 2 ]flt + 


c 2 - p \ 2 +| p* - p* | + \q-q\ 2 + \x dt 


(5.56) 


Where Ci and c 2 depend on the coefficients and the bounds on states and adjoints. Combining 
four of these estimates gives 
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^(P~ P) 2 (t f ) + p * - p*) 2 (t f ) + ^(q - q) 2 { 0 ) + (w-w) 2 \lt + ^(z - z) 2 ( 0) + 

h\\p - pf +| P* - P*f + \q~ +|z- z\ 2 \it < (c, + c 2 e utf ) 

| [p - pf +| P ~ P*\ +|y “ *f + |w - w| 2 + \z - zf] (5.57) 

0 

Thus from above equation we conclude that; 

(q + C 2 e Mf )j|p -Jjf +\p* - p r \~ +\q-q\ 2 + vv — vv| 2 + |z — Z " ]f// <0 (5.58) 

o 

Where C ] , C 2 depend on the coefficients and the bounds on P, P*,w,z 

~ ~ l ( A — C ^ 

If we choose A such that A>C l + C 9 and t f < — In\ — = — - , then p = p 

3A [ C 2 J 

p = p ,w = w;z = Z. Hence the solution is unique for small time 


6.1 Adomian Decomposition Technique (ADM) 

An explicit construction of approximate non-perturbative solutions of the system (5. 5-5. 8) was examining using Adomian 
decomposition method .The equivalent canonical form of this system is as follows: 

Tit) = 7(0) + £ sdt + (r - or)£ Tdt - JV + £ TT*dt - (1 - u x )k (1 + «,)"'£ V<#.£ VTdt + 


p\j*dt (6.1) 

7*(f) =T\0) + il-u l )k(l+ ay^Vdt.^VTdt -ij3 + p)\*T*dt (6.2) 

V(t) = V (0) + (1 - m 2 )(V/?£ T*dt - /£ Vdt (6.3) 

W{t) = W(0) + u 2 Np\ X Q T*dt - r [ Wdt (6.4) 

Adomian decomposition method was used to obtain the solutions of equations (3.290-3.293) are considered being as the sum of 
the following series; 

T = f J T n (6.5) 

n = 0 

r*|x (6.6) 

V = f j v n (6.7) 

n= 0 

W = Y j W n ( 6 . 8 ) 

n= 0 

Then we approximate the non-linear terms in the system as follows; 
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TT = Y J (\(T 0 ....T n , 

n = 0 

To-.-Tj) 

(6.9) 

VT = Y j (B n (V 0 ....V n , 

n = 0 

T 0 ....T n )) 

(6.10) 

7T* = X(c„(T 0 .... T n , 

T'~X)) 

(6.11) 




Iv 




( 6 . 12 ) 


(6.13) 


(6.14) 


The non-linear function A n , B n , C n , are called Adomian’s polynomials. Substituting equation 
(6.5-5.14) into (6.1-6.4) then we have; 


±T n = 7X0) + M + (r - «)£' fj„dt- - J- £ X- \dt + ~ 


fc(l - MlX 1 + XCE V n^[E B , dt + p[ Yjldt 

71=0 71=0 77=0 

(6.15) 

fx = r'(0) + *(i -«,xi + [Jy u d‘\ x „Y^. d > - (4+ 

(6.16) 

|)V„=V(0)+(1- u 2 )Np\lfj;dt - rl’jy.dt 

(6.17) 

Yw n = W( 0) + u 2 Nj3^^T*dt - 

(6.18) 

From equation (6.15-6.18) we define the following scheme: 
T 0 =T(O) + st 

(6.19) 

T* =T*( 0) 

(6.20) 

V 0 = V(0) 

(6.21) 

s 

II 

(6.22) 
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fr n = (r - a)\ X jy n dt - ^\X A n dt + j^!X C « dt - 

n=0 n = 0 1 max n=0 1 ^ n = 0 

+ ayiXvMXw + p\X T - dt 

n = 0 n=0 n=0 

(for n>0) 

IX = *(1 - Ml Xi + -(p+ p)\lfr:dt 

(for n> 0) 

=(i -u a )Np\lfr;dt- 7 \lfy n dt 

(for n> 0) 

fw n =u 2 Np\lfr:dt-y\lfw n dt 

(for n > 0) 

Using Equation (6.9-6.14). Some of the Adomian polynomials can be obtained as 

follows; 

F(t) = T 2 
We first set 

T = fr„ 

n=0 

Substituting (6.28) into (6.27)gives 
Fit) = (T 0 +T l+ T 2 +T 3 +T 4 +T 5 +...J 

Expanding the expression at the right-hand side gives 

F(t) = T 2 + 27’ 0 7J + 2T () T 2 + T 2 + 2 T 0 T 3 + 2f [T 2 + ... 

The expansion in equation (6.30) can be rearranged by grouping all terms with the sum of the 

subscripts of the components of the same. This means that we can rewrite equation (6.30) as 
Fit) = T 2 + 2T 0 f + 2 T 0 T 2 + T 2 + 2T 0 T 3 + 2 T 0 T 2 + 2 T 0 T 4 + 2fT 2 + T 2 + 

2 r 0 r 5 + 2fT 4 + it 2 t 3 + .... 

This gives Adomian polynomials for Equation (3.316) by 

A = T 0 2 
A = 277, 

^ = 2 T 0 T 2 +T 2 
A 3 = 2T 0 T 3 + 2T 0 T 2 
A 4 = 2T 0 T 4 + 2TJ 2 + T 2 
A 5 = 2 T 0 T 5 + 2 T x T 4 + 2 T 2 T 3 

Fit ) = VT 
We first set 


(6.23) 

(6.24) 

(6.25) 

(6.26) 

(6.27) 

(6.28) 

(6.29) 

(6.30) 

(6.31) 

(6.32) 

(6.33) 
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v = £v„ 

n = 0 

T = fj. 

n = 0 

Substituting (6.34 and 6.35) into Equation (6.33) yields 

F(T) = (v 0 +y +y 2 +y 3 +y 4 + v 5 +....)x(r 0 +t x +t 2 +t 3 +t 4 +r 5 +....) 

Multiplying the two factors gives 

m = V 0 T 0 + Ty + vy + T 0 V 2 + V 0 T 2 + Ty + T 2 V 0 + T 0 V 3 + 7 \v 2 + t 2 v 2 + t 3 v 0 + t 0 v 4 
+v 0 t 4 + t x v 3 + vy 3 + v 2 t 2 + ... 

Collecting all terms with the same sum of subscripts of the component T n , we can rewritten 
Equation (6.37) in the form 

Fit) = Ty + vy + t 0 v 2 + ry + t 2 v 0 + t 0 v 3 + ry 2 + ry + t 3 v 0 + 
t 0 v 4 + t 3 v 3 + t 2 v 2 + t 3 v ; +t 4 v 0 + .... 

Consequently, the Adomian polynomials are given by 

B 0 =T 0 V 0 

B i = ry +v 0 t 1 

B 2 =T 0 V 2 +Ty+T 2 V 0 

b 3 = t 0 v 3 + ry 2 + Ty + t 3 v 0 
b 4 = t 0 v 4 + ry 3 + t 2 v 2 + ry, + t 4 v 0 


Tit) = TT* 

We first set 

r = Zr„ 

71=0 

T' = lx 

77=0 

Substituting (6.42 and 6.41) into Equation (6.40) yields 

F(T) = (T 0 +t 3 +t 2 +T 3 +T 4 +T 5 +....)x(r*o +T\ +T\ +t* 3 +T\ +T 5 +....) 

Multiplying the two factors gives 

F(T) = T 0 T* + T*T | + T 0 T* + T*T 2 + T*T X + T*T 0 + T*T 3 + T*T 2 + T 2 T X + T*T 0 + T*T 4 
+T 0 T* + T*T 3 + Ty* + T 2 T* +.... 

Colleting all terms with the same sum of subscripts of the components T n , we can rewritten 
equation (6.44) in the form 

F(T) = T 0 T* + T*T X + T 0 T* + T*T 2 + T*T X + T*T 0 + T*T 3 + T*T 2 + T*T X + T*T 0 + 

T*T 4 + T*T 3 + T*T 2 + T*T X + T*T 0 + 

Consequently, the Adomian polynomials, are given by 

C 0 = T () T 0 

c 1 =t*t 1 +tX 

C 2 =T*T 2 +T*T 1 + T 2 T o 

C 3 = T*T 3 + T*T 2 + T*T X + T*T 0 

C 4 = T*T 4 + T*T 3 + T*T 2 + T 3 X + T*T 0 


( 6 . 34 ) 

( 6 . 35 ) 

( 6 . 36 ) 

( 6 . 37 ) 

( 6 . 38 ) 

( 6 . 39 ) 

( 6 . 40 ) 

( 6 . 41 ) 

( 6 . 42 ) 

( 6 . 43 ) 

( 6 . 44 ) 

( 6 . 45 ) 

( 6 . 46 ) 
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IV. Numerical results and discussion 



Fig 1: Uninfected CD4 + T-cell count with various a x and without treatment against time 
when s = 15;/? = 0.000024; ju x = 0.02; a = 500; c = 2.4; p 3 = 0.26; q = 0.1; 
k = 0.01; p = 0.05; 5 = 0.005; y = 0.01; p 2 = 0.40; a = 0.89 



Fig 2: Infected CD4 + T-cell count with various a ] and without treatment against time 
when s = \5\P = 0.000024; //, = 0.02; <7 = 500; c = 2.4; p, = 0.26; q = 0.1; 
n = 0.01; p = 0.05; 5 = 0.005; y = 0.01; p 2 = 0.40; « = 0.89 
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Fig 3: Viral particles with various a x and without treatment against time 
when s = \5\P = 0.000024; //, = 0.02; <r = 500; c = 2.4; ju, = 0.26; q = 0.1; 
n = 0.01; p = 0.05; 8 = 0.005; y = 0.01; p 2 = 0.40; a = 0.89 



time in (f-10 1 ) 

Fig 4: Uninfected CD4 + T-cell count with various CC X and treatment against time 


when s = 15\J3 = 0.000024; //, = 0.02; <7 = 500; c = 2.4; p, = 0.26; q = 0.1; 
n = 0.01; p = 0.05; 5 = 0.005; y = 0.01; p 2 = 0.40; « = 0.89 
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Fig 5: Infected CD4 + T-cell count with various <Z l and treatment against time 
when s = 15;/? = 0.000024; ju \ = 0.02; a = 500; c = 2.4; p 3 = 0.26; q = 0.1; 
n = 0.01; p = 0.05; 8 = 0.005; y = 0.01; p 2 = 0.40; a = 0.89 



time in wcvtfAs (MO 1 ) 

Fig 6: Infectious viral particle with various a x and treatment against time 
when s = 15] J3 = 0.000024; p x = 0.02; cr = 500; c = 2.4; p 3 = 0.26; q = 0.1; 

K = 0.01; p = 0.05; 8 = 0.005; y = 0.01; // 2 = 0.40; a = 0.89 
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Fig 7: Non- infectious viral particle with various a x and treatment against time 
when s = 15] J3 = 0.000024; //, = 0.02; cr = 500; c = 2.4; ju 3 = 0.26; q = 0.1; 
n = 0.01; p = 0.05; 5 = 0.005; y = 0.01; p 2 = 0.40; a = 0.89 


V. Conclutions 

In this paper, the study of the behavior of the mathematical 
model adopted to describe the HIV infection provides an 
explanation for the viral sudden deviation during the 
suppression of HIV. Adomian Decomposition Method 


(ADM) was also employed to compute an approximation to 
the solution of the non-linear system of differential equations 
governing the problem. Maple is used to carry out the 
computations and graphical results were presented to 
illustrate the solution. 


Appendix 

### ADOMIAN DECOMPOSITION METHOD ### 

>restart:a:=array (0. .50) :x:=array(0. .50) :y:=array(0. .50) : APROX : =array ( 0 . .50) :y[0] := 
0 . 0 : x [0] : =0 . 0 : h : =0 . 1 :M: =8 :N : =10 : a [0] : =719 . 57 :m: =0 :beta : =0 . 3 : alpha : =0 . 02 : r : =3 . 0 : g: =2 
. 4 :k:=0. 000027 :T [max] : =1500 . 0 : rho : =0 . 01 : s : =15 :b [0] :=27:c[0] :=3341:d[0] :=0:u[l] :=0.0 
083 :u [2] : =0 . 0083 : U0 : =s*t :U1 : =0 . 0 :U2 : =0 . 0 :U3 : =0 . 0 : TO : = (a [m] +t*dif f (a [m] ,t) )+(U0) :I0: 
= (b [m] +t*dif f (b [m] ,t) )+(Ul) : vO := (c [m] +t*diff (c [m] ,t) )+(U2) :w0 := (d[m] +t*diff (d[m] ,t) 
)+(U3) : A0 : =T0 A 2 : B0 : =v0*T0 : CO : =T0*I0 : T1 : =alpha*int (TO , t=x [m] . . t) + (r*int (TO , t=x [m] . .t 
) ) ( (r/T [max] ) *int (A0 , t=x [m] . . t) ) ( (r/T [max] ) *int (CO , t=x [m] . . t) ) ( (lu [1] ) *k*int (B0 , t=x 
[m] . .t) )+(rho*int(I0,t=x[m] . . t) ) : II := ( (1- u[l])* 

k*int (B0 , t= [x] . . t) ) ( (beta+rho) *int (10 , t=x [m] . . t) ) : vl : = (lu [2] ) *N*beta*int (10 , t=x [m] . 

. t) g*int (vO , t=x [m] . . t) : wl :=u [2] *N*beta*int (10 , t=x [m] . . t) g*int (wO , t=x [m] . . t) : A1 :=2*T 
0*Tl:Bl:=v0*Tl+vl*T0:Cl:=T0*Il+Tl*I0:T2:=alpha*int(Tl,t=x[m] . . t) + (r*int (T1 , t=x [m] . . 
t) ) - ( (r/T [max] ) *int (A1 , t=x [m] . . t) ) - ( (r/T [max] ) *int (Cl , t=0 . . t) ) - ( (1-u [1] ) *k*int (B1 , t 
=x[m] . .t) )+(rho*int(Il,t=x[m] . . t) ) : 12 : = ( (1-u [1] ) *k*int (B1 , t=x [m] . . t) )-( (beta+rho) *i 
nt (II , t=0 . . t) ) : v2 := (1-u [2] ) *N*beta*int (II , t=x [m] . . t) g*int (vl , t=x [m] . . t) : w2 :=u [2] *N* 
beta*int (II , t=x [m] . . t) g*int (wl , t=x [m] . .t) : A2 :=2*T0*T2+T1 A 2 :B2 :=v0*T2+vl*Tl+v2*T0 :C2 
:=T0*I2+Tl*Il+T2*I0:T3:=alpha*int(T2,t=x[m] . . t) + (r*int (T2 , t=x [m] . .t)) ( (r/T [max] ) *in 
t (A2 , t=x [m] . . t) ) ( (r/T [max] ) *int (C2 , t=0 . . t) ) ( (lu [1] ) *k*int (B2 , t=x [m] . . t) ) + (rho*int (I 
2 , t=x [m] . .t)) :I3:=((lu[l])*k*int(B2,t=x[m] . .t)) ( (beta+rho) *int (12 , t=0 . . t) ) :v3:=(l-u 
[2] ) *N*beta*int (12 , t=x [m] . . t) g* 

int (w2 , t=x [m] . . t) : w3 :=u [2] *N*beta*int (12 , t=x [m] . . t) -g*int (v2 , t=x [m] . . t) : A3 :=2*T0*T3 
+2*T1*T2 : B3 : =v0*T3+vl*T2+v2*Tl+v3*T0 : C3 : =T0*I3+T1*I2+T2*I1+T3*I0 : T4 : =-alpha*int (T3 , 
t=x [m] . . t) + 

(r*int (T3 , t=x [m] . . t) ) - ( (r/T [max] ) *int (A3 , t=x [m] . . t) ) - ( (r/T [max] ) *int (C3 , t=0 . . t) ) - ( ( 
1-u [1] ) *k*int (B3 , t=x [m] . . t) ) + 

(rho*int (13 , t=x [m] . . t) ) : 14 : = ( (1-u [1] ) *k*int (B3 , t=x [m] . . t) ) - ( (beta+rho) *int (13 , t=0 . . 
t) ) : v4 := (1- u [2] ) *N*beta* 
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int (13 , t=x [m] . . t) - g*int (v3 , t=x [m] . . t) :w4:=u[2]* 

N*beta*int (13 , t=x [m] . . t) -g*int (w3 , t=x [m] . . t) : 

Eq : =T0+T1+T2+T3+T4 : Eq2 : =10+11+12+13+14 : Eq3 : =v0+vl+v2+v3+v4 : Eq4 : =w0+wl+w2+w3+w4 : 

T = 719.57 + 1030.113124t + 34.97820628t 2 -719.6610060t 3 -41.31829380t 4 + 


T = 27 + 56.00465683t- 39.08 106622t 2 - 6.185455807t 3 +11.57469526t 4 +. 


V = 3341 -7937.173347t + 9608.980705t 2 -7735.947760t 3 +4631.970360t 4 +. 


w = .6733470426t-.1097224390t 2 -,2370917650t 3 +.1037093337t 4 + 
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